Early reionization by decaying particles and cosmic microwave background radiation 
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We study the reionization scenario in which ionizing UV photons emitted from decaying particle, 
in addition to usual contributions from stars and quasars, ionize the universe. It is found that the 
scenario is consistent with both the first year data of the Wilkinson Microwave Anisotropy Probe 
and the fact that the universe is not fully ionized until z ~ 6 as observed by Sloan Digital Sky Survey. 
Likelihood analysis revealed that rather broad parameter space can be chosen. This scenario will be 
discriminated by future observations, especially by the EE polarization power spectrum of cosmic 
microwave background radiation. 



I. INTRODUCTION 

The dark age between recombination and today has 
begun to be revealed to be very complicated by recent 
observations. The first witness is the Gunn-Perterson 
trough observed by Sloan Digital Sky Survey (SDSS), 
which implies that the universe was not fully reionized 
until 2^60. The second shed light on the begin- 
ning of reionization history, which may have started as 
early as z ~ 20, by the spectrum of cosmic microwave 
background (CMB) radiation observed by Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 2]. 

Conventionally the reionization is thought to be caused 
by the UV photons emitted from the early forming stars, 
where the ionization fraction of the universe is approx- 
imated by a step function ;33. In fact, the early reion- 
ization may be possible with large efficiency for star for- 
mation at high redshifts (z ~ 20) HH0. (See also 
for difficulty in early star formation with running spec- 
tral index suggested by WMAP.) However, it makes the 
ionization fraction unity well before z ~ 6 if it should ex- 
plain large optical depth Top ~ 0.17 suggested by WMAP. 
On the contrary, the Compton optical depth for reion- 
ization is not large enough if one considers UV photons 
from forming stars to be consistent with the data by 
SDSS 0,0. 

In order to avoid such difficulties, one should think 
more complicated reionization scenarios. For example, 
with the use of unconventional initial mass function and 
different photon emission processes between Population 
II and III stars, one can obtain the complicated ionization 
history which has a period of partial reionization before 
z ~ 6. In particular, some authors |^ 0, showed 
that the reionization occurs twice. 

Here we focus on another possibility. The idea is that 
the additional UV photons from decaying particles can 
explain the earlier reionization, while contributions from 
conventional stars and quas ars make the ionization frac- 
tion unity at z 6 ll3L . Although the existence 
of such particles may be suspicious, the ionization his- 
tory is very realistic once we admit it. The photons from 
the particle decay are characterized by three parameters: 
mass, abundance, and lifetime. In this paper, we concen- 



trate on the case that the mass is 30 eVt , which emits 
two 15 eV photons to ionize the universe. For the dif- 
ference between larger mass and this case, see We 
will see the decaying particle scenario is consistent with 
WMAP and SDSS observations. 

The structure of this paper is as follows. In the next 
section, we review how we obtain ionization history due 
to photons emitted from stars, quasars, and decaying par- 
ticles. Temperature (TT) and temperature-polarization 
(TE) power spectra of CMB observed by WMAP are used 
to evaluate the significance of this scenario in sect. HI. 
Section IV is for our conclusions. 



II. IONIZATION HISTORY 

In our scenario, the sources for ionizing UV photons are 
decaying particles, quasars and stars. Stars (and quasars) 
are indeed responsible for the full reionization at z ^ 6. 
On the other hand, contribution of UV photons from the 
decaying particles keeps ionization fraction much smaller 
before z ~ 6 for a long period. 

We follow the thermal history from z > 10'^ includ- 
ing recombination epoch, calculating the ionization frac- 
tions of hydrogen (HE) and helium (HeE and Hell), 
and the electron temperature, based on the argument of 
Fukugita and Kawasaki 15], where the hierarchical clus- 
tering scheme of the cold dark matter scenario is used. In 
addition, we include the sources of ionizing UV photons 
from decaying particles 12|. Some different treatments 
of stars and QSOs from 15] can also be found in [T^ . 

We assume that the particle emits two photons with 
monochromatic energy of half mass of that particle, i.e., 
E~f = The number density of (/)-particle is written 

as 

= n0(O)(H-z)3e"^, (1) 



t Here we consider this decaying particle as a part of cold dark 
matter (CDM). It can be regarded as CDM in such a way that 
a scalar condensate is coherently oscillating. 
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FIG. 1: Evolution of ionization fraction, = 10 sec (red), 
10^^ sec (green), and 10^'' sec (blue) are plotted. For each 
lifetime, both Top = 0.1 (solid) and 0.4 (dashed) are shown. 



FIG. 2: Likelihood contour of the decaying particle scenario. 



where is the lifetime of i/i-particle. Notice that those 
emitted photons with Ej > 13.6 eV can ionize hydrogen 
atoms. Then one can write a source term for the decaying 
particle as 



~dt 



T4, 



(2) 



In order to calculate how many photons are emitted, the 
abundance, mass, and the lifetime of the particle should 
be fixed. As mentioned above, we choose — 15 eV, 
since this case represents important features of decaying 
particle scenario. We calculate ionization histories with 
T0 = 10^'' — 10^^ sec (and 10^^ sec for some cases), adjust- 
ing the abundance Jl^ to get the desired value of optical 
depth Top, whose range we explored is Top = 0.1 — 0.4. It 
is defined by 
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where dr is the Thomson cross section and np is the 
electron number density for standard recombination. We 
subtract this term in order to estimate only the effect of 
reionization. 

Typical ionizing histories are shown in Fig. ^ Here 
ionization fraction is defined as Xe = n^/riH, where rig 
and uh are the number densities of electron and hydro- 
gen atom, respectively. 

As we will see later, cosmological parameters such as 
Hubble parameter, i?o, matter density, fim, baryon den- 
sity, rib, and the spectral index of primordial power spec- 
trum, ris, are varied such that they give the minimum 
for fitting both TT and TE spectrum observed by 
WMAP in the flat ACDM cosmology. 

As we can see in Fig. ^ ionization due to decaying 
particles endures for a long period, and the ionization 



fraction remains very small (10 



-3 < 



Xe 



< 



10-1) until 



z ~ 6. Even in these cases, the optical depth becomes 
very large. 



III. COSMIC MICROWAVE BACKGROUND 

Now we look for how these ionization histories are con- 
sistent with CMB observation by WMAP. To this end, we 
calculate the power spectrum with obtained evolution of 
the ionization fraction and matter temperature as inputs 
to the code modified from CMBFAST 1^, and evaluate 

using the code provided by WMAP |12| ■ We search for 
the range 10^^ — 10^^ sec, for the lifetime of the decaying 
particle, and Top = 0.1 — 0.4. For fixed lifetime and op- 
tical depth, the minimum is estimated adjusting flh, 
flm, EIq, Us, and the amplitude of the spectrum. As an 
indicator for the likelihood, we use Ax^ = x^ ~X^\hest 
make the contour of likelihood with grid-based analysis. 
Here we take X^j^g^j = 1428.6 as the six-parameter 
X^ minimum for ACDM model with WMAP-only data. 

The likelihood contours are shown in Fig. [5] (See also 
the Table HI below.) Here, we converted Ax^ into the 
confidence level for two degrees of freedom. One can see 
the ionization history with Top ^ 0.2 and = 10^'* — 10^^ 
sec is consistent with WMAP data. In general, there is 
little dependence on the lifetime of decaying particle. For 
larger optical depth, ^ VS^^ sec seems more preferable 
a bit, although its probability is very small because x^ = 
1434.7 and 1441.2 for Top ~ 0.3 and 0.4, respectively. We 
omit the lifetime longer than 10^^ sec from the figure, 
because the ionization fraction becomes unity before z ~ 
6 for Top ^ 0.3 even if there is no contributions of stars 
and quasars. Of course, if Top 0.2, it is consistent 
with observations as well. For example, x^ — 1432.6 and 
1430.5 for Top — 0.2 and 0.1, respectively. Notice that, 
for fixed optical depth, ionization histories with lifetime 
longer than 10^^ sec are identical to 10^^ sec case, since 
it is the value of /t^ that concerns with the amount of 
emitted photons; i.e., the abundance should be ten times 
larger for the ten times longer lifetime. This relation 
holds until the abundance reaches fi^ ~ r2„i. 

We should mention that x^ itself is smaller if reion- 
ization takes place instantaneously for the same optical 
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TABLE L Cosmological parameters which make minimum 
for typical optical depth and hfetime of decaying particle. We 
also show the contributions to the total ^ from TT and 
TE spectra. Lifetime is measured in unit of sec. Baryon 
abundance, fit., present Hubble parameter, Hq, and matter 
abundance, f2m are shown the ratios to best fit values for 
WMAP-only data: Vl^ = 0.047, i?o = 72 km/sec/Mpc, and 
= 0.29. 







TT TE 


Tls fij, Q,m 


10 0.417 
10" 0.316 

lU U.ZUo 

10" 0.098 


1450.1 
1438.9 

1429.8 


978.9 471.2 
976.7 462.1 

974.7 455.1 


1.2U 1.11 1.U2 U.99 
1.13 1.07 1.01 0.99 
i.uo i.uo u.yy u.yo 
0.99 1.00 0.98 1.00 


10^^ 0.405 
lO^'^ 0.314 
10^^ 0.202 
10^^ 0.097 


1438.8 
1434.7 
1431.5 
1430.4 


979.7 459.1 

978.8 455.9 
976.5 455.0 
975.5 454.9 


1.12 1.08 1.01 1.00 
1.07 1.06 1.00 0.99 
1.02 1.04 0.98 1.02 
0.97 1.01 0.96 1.03 


10^^ 0.407 
10^^ 0.307 
10^^ 0.201 
10^® 0.098 


1443.7 
1437.0 
1432.4 
1431.4 


986.0 457.7 
981.6 455.4 

978.1 454.3 
975.8 454.9 


1.05 1.05 0.99 1.00 
1.02 1.02 0.99 0.97 
0.99 1.02 0.97 1.01 
0.98 1.00 0.98 1.00 


10" 0.199 
10" 0.099 


1432.6 
1430.5 


978.6 454.1 
976.0 454.6 


0.98 1.03 0.96 1.03 
0.96 1.02 0.95 1.05 



depth. As we have mentioned several times, however, 
such instantaneous reionization histories do not respect 
the Gunn-Peterson trough at z 6 observed by SDSS. 
On the other hand, optical depth becomes ~ 0.04 if the 
ionization fraction at z 6 is suppressed enough to meet 
SDSS data. (When using our code without contribution 
of decaying particle, we get Top — 0.05.) 

Now let us see the CMB spectra in detail. We show 
TT, TE, and EE power spectra of CMB in Figs. EHHl 
For TT and TE spectra, we also plot the WMAP data. 
One can see from Fig. O that all the models look almost 
identical for £ > 100. The biggest differences can be seen 
in low multipole region. This happens because large op- 
tical depth reduces height of the acoustic peaks, which 
results in larger power in low multipoles (I ^ 100) when 
the height of the first acoustic peak is fixed. In addi- 
tion to this, the spectrum should be more tilted to the 
blue side for larger optical depth in order to fit the data, 
i.e., the well-known Top — degeneracy, which lowers 
the plateau. Thus the outputs are made by competition 
between the two opposite effects. 

Important point at present is that TE power spectrum 
does not have much ability of determining ionization his- 
tories using the first year WMAP data. The errors are 
still large, and even though those histories which do not 
seem to fit the TE data well have small (see Fig.0]and 
Table nj if the fit to TT data is good enough. This may 
lead to the fact that there is little dependence of the x^ 
on the lifetime especially for small optical depth region. 
Reversely, we could distinguish among different histories 
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FIG. 3: TT power spectra for typical ionization histories: 
= lO^** sec and Top ~ 0.1 (red), = lO'^^ sec and Top ~ 0.4 
(green), — 10^® sec and Top ~ 0.3 (blue), = lO^'^ sec and 
Top ~ 0.2 (purple), and instantaneous reionization which leads 
to Top ~ 0.17 (dashed black). In addition, WMAP data is also 
shown. 
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FIG. 4: TE power spectra for typical ionization histories for 
= 10^* sec (top panel), 10^® sec (middle panel), and 10^^ 
sec (bottom panel). In each panel, the lines are for Top ~ 0.4 
(red), 0.3 (green), 0.2 (blue), and 0.1 (purple). The instanta- 



neous reionization with To 



0.17 is shown in dashed black 



line, and the data is from WMAP. 



with various optical depth so long as the errors diminish. 

Although one may discriminate among the models us- 
ing TE spectrum, EE spectrum seems the most promising 
tool for that purpose. Most different feature of decaying 
particle scenario is the larger power at £ ~(a few) x 10 due 
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FIG. 5: EE power spectra for typical ionization histories for 
= 10^'* sec (top panel), 10^^ sec (middle panel), and 10^^ 
sec (bottom panel). In each panel, the lines are for Top — 0.4 
(red), 0.3 (green), 0.2 (blue), and 0.1 (purple). The instanta- 
neous reionization with Top — 0.17 is shown in dashed black 
line. 

to enduring UV photon emission. Moreover, we may dis- 
tinguish even the lifetime of decaying particle depending 
on whether there is a dip around £ ^ 10 and the steep- 
ness of the slope below i ~ 100. Therefore, EE spectrum 
expected to be seen by Planck satellite or even WMAP 
should (dis)prove the photons from decaying particles in 
the near future. 



IV. CONCLUSION 

We have sought for a natural reionization history which 
respects both WMAP and SDSS observations. Decay- 
ing particles provide UV photons which reionize the uni- 



verse from rather early period as suggested by WMAP, 
but keeps the ionization fraction very small until usually 
considered UV sources such as stars and quasars bring 
up the fraction to unity abruptly at z 6 as observed 
by SDSS. The central concern is how well such scenario 
goes on. To this end, we have calculated TT and TE 
power spectra of CMB, and done the likelihood analysis. 
We have found decaying particle scenario is consistent 
with both observations in rather broad parameter space. 
Even particles with lifetime lO^'* sec, corresponding to 
the redshift z ~ 280, which have the abundance to give 
Top ^ 0.3, seems good as well. 

Looking into the detail, it is the TT power spectrum 
that determines the value of dominantly. This is the 
reason why the reionization history, whose TE spectrum 
does not seem to fit the data so well, can have small x^j 
such as 1429.8 for Top ~ 0.1 and = 10^* sec, for ex- 
ample. When the quality of TE data improves, then one 
may be able to tell which reionization history is prefer- 
able. As for the ability of discriminating each history, 
EE spectrum may be the best tool. Decaying particle 
scenario has a peculiar feature that there is much more 
excess of the power in the range 20 ^ £ ^ 100 compared 
with instantaneous reionization case. Moreover, it may 
be possible to observe even the lifetime of the decaying 
particle. 

Finally, we comment on a maybe serious problem if one 
regards early forming stars as the whole source of ioniza- 
tion photons because of the less powers on smaller scales 
observed by, say, 2dfGRS and SDSS. Since the decaying 
particle scenario does not owe to early forming stars as 
the UV source, it is completely free from such a problem. 
We have checked that our model works also for the run- 
ning spectral index suggested by WMAP. Since all the 
astrophysical UV sources may suffer this problem, the 
decaying particle could be the better source of ionizing 
UV photons. 
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